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1. Introduction 

Factorization underlies any theoretical prediction at hadron colliders, since it allows to sep- 
arate the short distance, perturbatively calculable physics from the non-perturbative ingre- 
dients such as the parton distribution functions. First arguments in favor of factorization 
for hard QCD processes to all order in perturbation theory appeared over three decades ago 
H|, H|, [|. Power counting, pinch analysis and physical pictures have been applied to 
analyze generic loop integrals at high energies in Refs. j|, f?J, which led to the seminal work of 
Collins, Soper and Sterman on proofs of factorization theorems to all orders in perturbation 



series (see Refs. ||, ||, 10, 11] and references therein), and has been accomplished for many 
processes of interest. In these proofs the Landau equations fi"2fl , and their physical interpreta- 
tion using the Coleman-Norton Theorem [O] , were used to identify the infrared singularities 
giving rise to the long distance physics in arbitrary Feynman diagrams. 

Singularities in Feynman diagrams arise if the integral over the loop momenta leads to 
pinched singularities, where the contour of the loop integration cannot be deformed to avoid 
the singularities in the integrands. The Coleman- Norton Theorem allows to map these pinched 
singularities onto configurations containing on-shell particles. For most cases, pinched singu- 
larities are due to internal propagators becoming ultrasoft or collinear to external particles, 
but there are also pinched singularities associated with so-called "Glauber" 0, ^] modes 1 . 
These are internal modes which have transverse momentum much in excess of their longitudi- 
nal momentum, and can therefore not be thought of as being on the mass-shell. In particular, 



1 In Ref. M the same momentum region was called "Coulomb" 
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Glauber gluons connecting the remnants of the beams in hadronic collisions have been argued 
to threaten factorization [14, The presence of this transverse interaction causes signifi- 
cant complication in the proof of the factorization for the Drell-Yan process. However, after 
summing over the unobserved hadronic final states, the effects from Glauber gluons cancel, 
and factorization holds. This was explicitly shown to one loop order in and generalized 



to higher orders in Refs. [17, [l{J. Another approach was developed in Ref. [19], where the 
cancelation of Glauber gluons in the inclusive cross section was shown using the light-cone 
ordered perturbation theory. 

Recently, a different approach to factorization proofs has emerged 
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34, 35, 36, 37, 38], which is based on effective field theory 



techniques. Effective field theories describe the long distance physics using a limited set of 
degrees of freedom, with all short distance physics being integrated out of the theory and 
contained in the Wilson coefficients of operators. The relevant effective theory for collider 
physics is soft-collinear effective theory (SCET) (3^, 40, 41, , which in its usual formulation 
contains collinear and ultrasoft degrees of freedom. Factorization in SCET can be shown in 
a very straightforward fashion, since one can show at the level of the Lagrangian that the 
collinear and ultrasoft degrees of freedom decouple in SCET. The advantage of the SCET 
approach to factorization theorems is that all long-distance ingredients of the factorization 
theorem are defined as matrix elements of effective theory operators, which allows to use 
renormalization group equations to resum large logarithms that arise in most perturbative 
expressions. 

In its traditional formulation, SCET does not include Glauber gluons. Given their im- 
portance in traditional factorization proofs it is therefore crucial to understand if and how 
Glauber gluons enter the effective theory framework. An attempt to include Glauber gluons 
into SCET was made in 43], where the factorization of the DY cross section in the presence 
of a Glauber mode was reconsidered. However, there are several flaws in the arguments pre- 
sented by these authors, in particular their argument about the necessity for Glauber gluons 
in SCET fails to consider the overlap between ultrasoft, collinear and Glauber modes. As 
we will show in this work, a proper treatment of this overlap is crucial to understand the 
contribution of Glauber gluons. Note that the inclusion of Glauber gluons into the SCET 
Lagrangian has also been used to describe jet broadening in dense QCD matter |44|, [4q| . 

The main purpose of this paper is to investigate whether Glauber gluons are required 
in SCET or not. We consider a well known operator in SCET, the operator producing two 
energetic back-to-back quarks, and perform the matching calculation determining its short 
distance Wilson coefficient using two different choices for the external states. Given that 
the short distance physics has to be independent of the external states, a consistent effective 
theory has to give the same result for both of the calculations. We will show that SCET^j , the 
traditional formulation without Glauber modes, gives different results for different external 
states, while a theory that includes Glauber gluons, SCETq , will give the correct result 
for both choices. This unambiguously shows that Glauber gluons need to be included for a 
certain class of processes. It should be pointed out, however, that Glauber gluons may be 
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integrated out of the effective theory, since they can not be on their mass-shell and cannot 
appear as external particles in perturbation theory, leading to a potential between pairs of 
collinear fields in opposite directions |4q| . 

Having established the requirements of Glauber gluons in SCETg , we will study in some 
more detail the relationship between pinched surfaces and the effective theory. Using a simple 
graphical representation of the pinched surfaces one can easily understand the necessity of 
Glauber gluons in the matching calculation under consideration. We then proceed to use this 
analysis to derive what final states are necessary to require Glauber modes to give a non-zero 
contribution. 

The plan of the paper is as follows. In Section |2| we perform a one-loop matching calcu- 
lation within SCET^g involving DY amplitude topologies and show that the effective theory 
breaks down, while SCETg passes the consistency check. In Section || we perform the pinch 
analysis of a diagram with pure spectator interactions and identify the right modes of an 
effective theory for this process; then we discuss other processes where the Glauber mode 
plays a role. Finally we conclude in Section ^[ 



2. An explicit matching calculation 




Figure 1: One-loop examples of (qq\qTq\0) (left) and (^^IqTqlqq) with active-active, spectator- 
active, and spectator-spectator interactions. 



In this Section we will consider the well-studied SCET current O2 = XfJ^Xn- The Wilson 
coefficient is usually calculated by an explicit calculation using partonic external states of free 
back-to-back quarks with large energy 

(qq\qTq\0) = C2(qq\O2\0) + power corrections. (2-1) 

where O2 is built of collinear gauge invariant fields Xn = VFn£ n - The Wilson line W n contains 
n— collinear gluons and £ n is a two-component spinor describing an n— collinear quark field. 
(For more details see e.g. Refs [ p0| , p| , p2] , 23, 47] and References therein). The relevant 
Feynman diagram in the full theory is shown on the left of Fig. []]. 

It is a well known consistency requirement of effective theories that the short distance 
Wilson coefficient of any operators has to be independent of the long distance physics in the 
process, and in particular it has to be independent of the external states chosen for the match- 
ing calculation. Inspired by the fact that Glauber gluons are known to manifest themselves 
through interactions with beam remnants, we intentionally choose more complicated external 
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states, namely 

(Yl*\qrq\qq) = C 2 {^*\0 2 \qq) , (2.2) 

where the two photons in the initial state and the two quarks in the final state are back-to-back 
with large energy. In this matrix element, each of the two off-shell photons converts into a qq 
pair, with a quark and an antiquark from each photon are then annihilated by the operator, 
while the other two quarks end up in the final state. The required Feynman diagrams in the 
full theory are depicted in Figure [j]. We refer to the quarks that are annihilated by the operator 
as active, while the quarks that end up in the final state are called spectators, and divide 
the loop diagrams into active-active, active-spectator and spectator-spectator diagrams, as 
indicated in the Figure. Note that beside the internal propagators that are not involved in the 
loop integral, the active-active diagram is the same diagram as in the matching calculation 
with a qq in the initial state, and should by itself give the correct result for the matching 
calculation. This implies that the sum of the contributions to C 2 from active-spectator and 
spectator-spectator diagrams has to yield zero in a consistent effective theory. 

For the purposes of this paper, i.e. identifying contributing modes, it is sufficient to 
study the singularity structure of the integrands, and we will therefore omit the numerator 
structure of propagators for simplicity. In the rest of this section we proceed as follows: 
The Wilson coefficient C 2 is determined via a matching calculation of the simple matrix 
element given in Eq. fl2.ip . Then we repeat the procedure for the more complicated matrix 
element given in Eq. fl2.2| ). First, we will calculate the contributions to the DY process in the 
full theory straightforwardly. Next, we perform the matching calculation by computing the 
corresponding diagrams in SCETq , finding that C 2 is entirely reproduced by the active-active 
topology, as expected. The spectator-active diagrams in SCETq match the corresponding 
diagrams in the full theory. In the spectator-spectator topology, however, only the soft 
diagram in SCETq is of leading power, and this diagram does not reproduce the full theory 
diagram. It appears that this contribution would change both the Wilson coefficient and its 
anomalous dimension from the findings of Eq. ( |2.1| ), which is inconsistent. Finally we consider 
SCETq , which includes Glauber gluons. Repeating the matching once again, we find that 
Glauber gluons do not contribute for both Eq. ( |2.lD and the active-active part of Eq. (|2.2| ), so 
that C 2 remains unchanged. Interestingly, the active-spectator topologies in the full theory 
and SCETq again match identically due to the fact that the naive contribution from Glauber 
gluons is precisely cancelled by its overlap with the collinear gluon exchange. Finally the 
spectator-spectator topologies also match, and SCETq passes the consistency check. We 
emphasize that the full theory diagram is identically reproduced by the sum of the soft and 
the Glauber gluon exchange, and not by the soft gluon exchange alone. 

Before we consider the various diagrams in turn, let us set up some notation. At tree 
level we define 

{qqWq\Q)tree = (gg|O 2 |0)tree = 1 • (2.3) 
For the matching calculation denoted in Eq. ( |2.lD we denote the full theory matrix element 
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at one-loop by 

(qq\qTq\0) loop = I 3 , (2.4) 

where ^3 denotes the result of the 1-loop calculation. In conventional SCET^ , the gluon is 
either collinear or ultrasoft (which we for brevity call soft everywhere below). They scale like 
(A 2 ,1,A) and (A 2 , A 2 , A 2 ) in their light -cone components, respectively, where A is the small 
expansion parameter of SCET, and the light-cone components are as usual with respect to 
the null vectors n and n in the ±z direction. To one loop order the amplitude for this effective 
theory can therefore be written as the sum of tree level and three effective theory diagrams 

(^|O 2 |0) = l + /3 C + / 3 E + /f. (2.5) 

Thus, the Wilson coefficient is given by 

_ (qq\qrq\0) _ e 

2 "7^W 3 " ( 3 3 3) ' ( } 

For the alternative external states, Eq. (|2.2j), we assign momenta p + p and q + q to the 
two initial virtual photons, and p and q to the two outgoing(spectator) partons. The tree-level 
amplitude is simply given by the two propagators of the active quarks, such that we find in 
both the full and the effective theory 

(7*7*1 q^q \qq) tTee = (7*7*1 2 \qq) tree = ^2 • {2 ' 7 ) 

At one loop order, the full theory amplitude of the matrix element of our interest is equal to 
sum of three topologies, which we write as I3, I4, ^5, namely a triangle graph in the active- 
active, box graphs in the spectator-active and a pentagon graph in the spectator-spectator 
topology. We write for the full theory calculation 

(7*7*1 qTq \qq) = 4=2+ 4=2^ + 4 7 4° + 4^ + ^ > ( 2 - 8 ) 
pZqz pAgz qz pz 

where the prefactors of 1 /p 2 , 1 /q 2 take into account the propagators that are independent of 
the loop-momenta. The (n) and (ra) superscript on the I4 denote if the active quark is in the 



n or n direction. Note that the ^3 denotes exactly the same integral as in Eq. (2.4). 
In conventional SCET^j the one loop amplitude can be written as 

(7*7*1 2 \qq) = + ^(J| + jf + If) + ^(if* + + 8 + lf )s ) + If , (2.9) 

where the additional superscript c, c and s on the integrals denotes if the gluon is collinear 
in the n direction, collinear in the n direction or soft. Note that we have already used the 
fact that some graphs, i.e. j| , l| , If, If, are power-suppressed. The Wilson coefficient is 
therefore 

^2 



c 2 = i + i 3 -(ii + ii+m+p 2 



T (n) I An)c T {n)s 
J 4 \ 1 A 4 



2-2 



T {n) ( T {n)c T (n)s 
J 4 \ 1 4: "T J 4 



+p^[I 5 -II] . (2.10) 
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In theory SCETg , which adds Glauber gluons with momentum scaling (A 2 , A 2 , A) to the 
conventional SCET^j , there are new Feynman diagrams present at one loop. Denoting the 
triangle, box and pentagon diagrams with Glauber gluons in the loop by /f, I^ 9 , I^ 9 and 
Ii, respectively, we can write 



(7*7*1 2 \qq) 



G ~ 1-1 



- + + 4 + n + id + j 2 (it y + it )s + it )9 



+i(/f )e ' + if )s + if )a ) + ii + 1 9 



(2.11) 



such that we can write 



c 2 = i + i 3 -(ii+ii+ii + i 9 ) +P i 



r{n)c' An) 



+ ir + 1 



( n )g 



r(fi) 



(n)c' 



+ if )S + I 



( n )9 



+ pV[/ 5 -(/! + /; 



(2.12) 



Thus, the condition that the Wilson coefficient C 2 has to be independent of the external 
states chosen allows us to determine if the theory with or without Glauber gluons is correct. 
If SCET^j were correct one would find 



P 



r(n) 



r (n)c , T (n)s 

A T J-A 



+ r 



r(n) 



lf YC + I 



P 2 f[h 



(2.13) 



while the presence of Glauber gluons changes this condition to 



r(n)c' T {n) 



+ IT' + 1 



(n)g 



An) 



r (n)c' T (n)s 



+ IT' + I 



( n )a 



+P 2 ? Ih - (I! + if)] = . 



(2.14) 



We will explicitly show below that SCET^ does not satisfy its consistency check Eq. ( p. 13 ), 
while the effective theory which contains Glauber gluons does satisfy Eq. ( [2.14 ). We will 
see that this difference between SCET^ and SCETq happens only in the spectator-spectator 
topology, i.e. in the last terms of Eq. ( 2. 13| ) and Eq. (2.14). 

It is important to note that the collinear integrals differ between the two different versions 
of SCET. It is well known that care has to be taken when defining the contributions arising 
from collinear gluons, since the various modes in the effective theory have overlap regions 
which can lead to double counting. When integrating over collinear loop momenta in SCET^j , 
part of the integration is over a region in which the collinear momenta become soft. This would 
double count the soft region and therefore has to be removed from the collinear diagrams. This 
procedure is called "zero-bin subtraction" [ 48], and the collinear contributions in Eqs. ( [2.10 ) 
and (gig ) are all zero-bin subtracted. In particular we have 



(I, 



k)0s ) 



(I 



fc)0s , 



(2.15) 
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where i| denotes the naive unsubtracted integrals. When adding Glauber gluons, there are 
more overlapping regions, requiring a more involved zero bin subtraction. We can write 



1 k 



1 k 



k 



1 k 



{Ik )0ff0s 
(-^fc )0g0s 



)og + ( J k )os 
( J k )og + i J k )os 

II (Wo* ■ " (2-16) 

2.1 Full Theory one loop calculation 

The active-active topology giving rise to I 3 in the full theory is simply a standard scalar 
triangle integral. In d = 4 — 2e dimensions one finds 2 

d d l 



2,,2e 



h = H)<?V 
l 



r d a 
J W 



) d [Z 2 + iO] [{l+p) 2 +i0} [(l-qf + iO] 



4-7T p + q 



7T 2 p 2 

-5- +ln-— 

3 p+q 



In 



+ 0(e, A 2 ) 



The spectator-active topology can be written as the box integral 



(-%V £ 
1 



d d l 



{2ir) d [I 2 + iO] [{I - p) 2 + iO] [(I + p) 2 + iO] [(I - q) 2 + iO] 



(2.17) 
(2.18) 

(2.19) 



a 



Att q p 2 p + + p 2 p + 



^-2Li 2 
3 



2 -4- 

p2p+ 



+ 



In 



-2 -t 

p2^ 



— Z 7T 



7- 



In 



p-Tp^ _|_ p+p 2 

g 2 (p + p) 2 p+p 2 




+ O(e,A ). (2.20) 



An analogous expression is valid for the second spectator-active integral i^. 

The spectator-spectator topology giving rise to 1$ in the full theory can be calculated via 
a pentagon integral which by standard procedures can be reduced to sum of five box integrals. 
The result is 3 

* d l 1 



h 



2,.2e 



/• d a 

y (2^ 



= — M + M~ 

47T 



+Z7T In 



)d [p + i0 ] [(/ - p) 2 + iO] [(/ + p) 2 + iO] [(/ - q) 2 + iO] [(/ + q) 2 + iO] 
1 



p+p+{p+p) 2 q q (q + q) 2 



i_ o 2 

m — : — K In 



p-Tp* 



q q' 



p 2 p 2 q 2 q 2 



p+p+q-q~(M+M-) 2 



+ 



2m 



p+p+(p + p) 2 (M-) 2 -q-q-{q + q) 2 {M+) 2 



q q {q + qf 



p+p+(p + p)< 



(2.21) 



In all diagrams in this paper we assumed that off-shellness is positive: p ,p ,q ,q , (p + p) , (g + q) > 0. 
3 This pentagon loop integral and also integrals Jf , , (/| n ' c )o g have been calculated for a simplified case 
of "_L"-less kinematics: px = px = qx = qx = 0. Since we are free to choose any external states in the 
matching calculation we can easily achieve these conditions. 



- 7 - 



where we have defined 



M + = p + + p + , 



M~ 



+ ■ 



(2.22) 



To determine the Wilson coefficient C2 we need to subtract the effective field theory 
diagrams from this result, which should cancel all IR divergences and reproduce the scalar 
one-loop result 

rv. T 1 In -f 

C 2 



47T 



+ in l 



A* 



+ - In' 
2 



+ m In 



A* 



7 



-7T 



(2.23) 



p-rg- 12 

2.2 Conventional SCET: soft and collinear gluon exchanges 

The loop integrals for the effective theory modes in each topology can be found trivially by 
expanding the integrands in the full theory with the appropriate scaling of the gluon momenta 
before integrating over the loop momentum. The zero-bin subtraction integrals can be found 
similarly by expanding the effective theory loop integrals with the scaling of the overlap mode. 
All zero-bin integrals are scaleless in this effective theory, so they just convert infrared poles 
into the ultraviolet ones. 

The results for all diagrams in SCET^ are given by 



1 3 



i 3 



J 3 



a. 



1 



An p+q 
1 



1 



a 



4n p + q 
1 



In ^ + in 



In 4y + in 



In — + in 



+ 



7T 

12 



In + in 



+ 



TT 

12 



a 



4n p + q 



11 n 2 p + q~ 

-9 + - In — ^zjt 
e e p z q L 



ITT . 

+ T + 2 ln 



1 . 2 ft 2 p+q- 



2 -1 



+ inln 



p 2 q 2 



7T 

T 



T {n)c C^s 

-M — ~. — 



, p 2p+ . . 

In + in 



7n 2 



An q p 2 p + + p 2 p^ 



-2Li 2 - 



p 2 p + 

p2p+ 



r (n)s 



2 2^ + )2fc+)2 p 2 p 

+«7T In ,„ , q 1 \ ri ri + hi 



(p 2 p + + p2p+Yp 2 



p2p-\ 



In 



(p 2 p + + p 2 p+) 7 

(p + p) 2 p+p+p 2 2 p 2 p 2 p + 



1 4 ^ 
In M P 



1 



/I = — 



1 



4-7T p 2 p + + p 2 £T 



In 



2f!+ 



^ l ln gV^ ln fJ- 4 P+P+(q ) 2 



— Z7T In 



a. 



M+M~ 



e 2 

H 2 p 2 (p+) 2 q~ 
q 2 (p 2 ) 2 p + 

2 in 



p2p-\ 



3 2 



47r p + p + (p + p) 2 q q (q + q) 2 

+in In 



+ In - , ^ „ In 



p+p2 



q q 

q-q 1 



p 2 p 2 q 2 q 2 
p + p + q~q~ fi 4 



+ 3vr 2 



(2.24) 
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We can combine the results from the full theory and SCET^j calculation to find 



I3-I3-I3- 1 3 



P A [h-l 



(n)s r(n)c 



Us 
47T 



I In 



p + q~ 1 1 1 2 H 1 i-i A 4 '2 

i — i 1 — in — - — l-nrm — — 7r 

e 2 p+g j>+g~ 12 



0, 



p 2 g 2 



+ 



47r I p + p + (p + p) 2 q~q~(q + <?) 2 
27rip 2 g 2 

p+p+(p + p) 2 (M-) 2 -q-q-{q + q) 2 {M+f 



2tti 9 , u 2 

2vr 2 + 2vr i In — r 

e M+M- 



ln 



(A/~) 3 M+ 



(M-) 2 ^!l g "~ ((?+g " )T 



q-q-{q + qf 



(Af 



ln (M+)3M- 
+ ^2 p+p+(p+p) 2 
p+p+(j) + p) 2 



(2.25) 



As discussed before, the triangle diagrams are the same as with the external state, such 



that the combination I-i — II 



J 3 



K reproduces the known Wilson coefficient. While the 



box diagrams are the same in the full and effective theory and therefore the combination 
I4 — I^ s — I^ c is equal to zero, this is not true for the combination I5 — i|. This implies 
that the Wilson coefficient C2 is different when calculated with the two off-shell photons in 
the initial state, which shows clearly that SCET<j does not reproduce the infrared physics of 
the full theory. On top of that, we can see that the combination I§ — J| is UV divergent as 
well, indicating another failure of SCET<^ . 



2.3 SCET including Glauber gluons 

Having found a serious problem in SCET^g we now repeat our matching calculation including 
Glauber gluons. One interesting aspect of this theory is there exists a zero-bin contribution 
in the spectator-active topology that does not vanish, such that it is more than merely a 
"pull-up" contribution [48[| . The non- vanishing overlap is 

i7T ■ 1 ( p + p 2 +p + p 2 . n \l , 00fiA 

Z7rln T -— z-T-^0) ■ (2.26) 

e V V- 2 {p + +P + ) J\ 



Kh k 



Og_ J- 1 

47T q~ p + p 2 + p + p 2 



Continuing with the calculation of the other contributions, we find that several of the 
integrals are either equivalent to the SCET^j case or zero 

Il'=I c 3 , II' = 11, 4 = 0. (2.27) 
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For the remaining integrals we find 

Jn)c' _ Jn)c /r(n)c\ 
1 i — 1 4 — l J 4 )0g > 



r<»9 



- (4 n)C )og > 



If 



47T 1 p+p+(p + p) 2 g g (<? + <?) 2 L 6 

2vri M+M 

+ 



2vr 2 + 2tt z In 



M+M- 



+ p) 2 (M-) 2 - q~q-(q + g) 2 (M+) 2 



(AT 



In iMZl^t 
q-q-(q + q) 2 



(M H 



ln (M+)3Af- • 

2 p+p+(p +p )2 

p+p+[p + p) 2 



(2.28) 



The corresponding contributions to the Wilson coefficient Ci from different topologies 
are equal to 



h-n -it 



II - if = ^ 

3 3 4vr 



1 

- + 

e z 



ln 



— h-ln t _ + «vr ln — 

e 2 p + q~ 



7 



7T 

p^g" 12 



p 2 (/ 4 -4 n)c '-4 n)ff -4 n)s ) =0, 

A 2 (/ 5 -/|-/f) = 0. 



(2.29) 



Thus, SCETq with the inclusion of Glauber gluons does reproduce the correct Wilson coeffi- 
cient, which can be taken as a strong indication that it is the correct effective theory. 



3. Pinch analysis and power counting 

3.1 The spectator-spectator interaction in Drell-Yan 

As we showed in the previous Section, the spectator-spectator contribution required a Glauber 
gluon for the effective theory to reproduce the full theory result. In this Section we will 
investigate the structure of this contribution in more detail in order to get more insight into 
the required modes. The scalar integral contributing to the spectator-spectator diagram is 
given by 



H) 



dH 



(2tt) 4 I 2 + iO (I + p) 2 +i0(l- p) 2 + i0{l + q) 2 + i0{l- q) 2 + iO 



(3.1) 



Decomposing the loop momentum I into its light-cone components we arrive at the following 
form, which is suitable for first integrating over the + component by contours and leaving 
the _L components as a final integration: 



T < .xl f d 2 l ± fdl- AT _ n . [dl+Jj 



i 



(3.2) 
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Figure 2: A few representative "constellations" of singularities during the integration over l±. We 
show the three non-trivial intervals in l~ (regions (II), (III), (IV) above) and one sample projection 
onto the complex l + plane for each region. The poles in the complex l + plane can be pinched in 
different locations, depending on the values of l~ and l±. For example, in region (III) l~ counts as 
0(A 2 ), and there are two pinched singularities: first, there exists a value of l± for which l + is pinched 
between 23 and 24 ~ 0(A 2 ); second, there exists a (generally different) value of l± for which l + is 
pinched between z\ and 22 ~ 0(1). 



where [N~(l~)] = +p~){l~ — p~){l~ + q~)(l~ — q~)- The positions of the singularities 
Zi in the complex l + plane are functions of l~ and l±, as well as the external momentum 
components. Explicitly they are given by 

zq{1 ,l±) = — p — , 

I + p I +q 

: z hp , Z4{1 ,l±) = ; = h<? • (3.3) 

I — p I — q 

Note that the locations of the poles above or below the real axis changes during the integration 
over l~ at the transitions l~ = -q~ ~ 0(1), q~ ~ 0(1), ~P~ ~ C(A 2 ), p~ ~ C(A 2 ) and 
Z~ = 0. For Z~ < — q~ or Z _ > q , all poles are either above or below the real axis, such 
that the total integral vanishes. We can divide the remaining range of l~ into three regions 
— q~ < l~ < —p~, —p~ < l~ < p~ and p~ < l~ < q~ , as depicted schematically in Figure ||. 
We begin by considering the region — p~ < l~ < p~, labeled (III) in the Figure. In this 



z\(l ,l±) = 
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case the poles in Z\ and Z3 are below the real axis, while the poles in z<i and 24 are above the 
real axis. The pole in zq is either above or below the real axis 4 , depending on the sign of l~. 

Since l~ ~ 0(X 2 ) the pole location Zi, and hence the magnitude of l + after its residue 
is taken, depends only on the size of l±. In particular, zq ~ 0(/^/A 2 ) depends strongly on 
the power assignment of l±, which demonstrates the shortcomings of the cartoon in Figure ^. 
We find that a new approach to illustrate the l± dependence is in order, which we present in 
Fig. g. 

The top left part of Fig. || shows the magnitude of the pole position Zi as a function of l± 
on a double logarithmic scale. The green, orange, red, magenta and blue curves correspond 
to zq, zi, Z2, Z3 and Z4, respectively. We also visualize the sign of the pole by showing the 
Zi above the real axis by solid lines, those with poles below the axis by dashed lines, and zq, 
whose pole that can change location, by the dot-dashed line. Whenever a solid and a dashed 
line are present at the same point, the two corresponding poles can be pinched. For example, 
for l± ~ 0(A) there is a pinch between 2:3 and Z4 leading to l + ~ A 2 , but also a pinch between 
zq, z\ and Z2 leading to / + ~ 1. 

Next, for a given l±, we can also determine the magnitude of the contribution of each of 
the residues to the integrand of 1$ . In the top right part of Fig. || we illustrate the magnitude 
of the integral, taking into account the scaling of the d A l. In this Figure we have chosen to 
close the contour above the real axis, so that the poles below the real axis do not contribute. 
Note that leading power corresponds to A -4 , consistent with Eq. (|2.8| ). 

These diagrams can be used to identify the modes required in an effective theory. From 
the top right plot we see that a leading order contribution A -4 can only come from the residue 
of 24, and that l± has to be between A and A 2 . The residue from zi always leads to a power 
suppressed contribution, in agreement with the explicit calculation performed in the previous 
Section. From the top left plot we see that the corresponding poles are indeed pinched and 
we also read off that for these cases the momentum component l + scales as A 2 . Thus, the two 
scalings of the loop momentum that give rise to leading order effects in the spectator-spectator 
contribution are 

P ~ (A 2 , A 2 , A 2 ), and P-(A 2 ,A 2 ,A), (3.4) 

which are the soft and Glauber mode. We will come back to the possibility of I 1 - being 
somewhere in between A and A 2 later. 

So far we have only considered the case l~ ~ 0(A 2 ), and in particular the region l~ € 
[— p~ ,p~]. Outside that region one of the poles z\, Z2 will cross the real axis onto the other 
half plane, unpinching a singularity, which only results in the removal of the already power- 
suppressed (collinear) contribution of its residue. For completeness we shall also mention that 
the case l~ ~ O(l) will again lead to no further leading-power contributions. 

It is an easy exercise to draw plots analogous to the top right plot of Figure || for the 
effective theory. As one obtains the integrals in the effective theory by expanding the full 
integral about the loop momentum modes, one similarly finds that the top right plot in 

4 In the Figure we chose to indicate this by placing the pole directly on the real axis. 
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Figure 3: Top Left: magnitude of pole locations as a function of l±. Dashed lines denote poles in the 
lower half plane, while solid ones are in the upper half plane. The dash-dotted line may be on either 
side. Top Right: magnitude of the residues of poles in the upper half plane (including integration 
measures). The color coding is identical to the one on the left. Lower Row: magnitude of the residues 
for the effective modes soft (left) and Glauber (right). 



Figure || is simply expanded around the modes, which is shown in the lower row of the 
Figure. In this case, expanding the integrand around soft loop momentum will reproduce the 
plot below l± ~ 0(A 2 ) and continue the straight segment to the right of it indefinitely. This 
is shown in the lower left plot. Thus, the soft contribution reproduces the infrared region 
and picks up a 1/euv singularity. On the other hand, an expansion around the Glauber 
scaling (lower right plot) will reproduce the region around l± ~ A and pick up a 1/eir pole. 
The overlap between those two modes is obtained by expanding the Glauber integral around 
the soft limit, which is a scaleless integral and provides the subtraction of both poles, thus 
reproducing the full theory at leading power. 

In the above we have only considered the two discrete choices l± = A or lj_ = A 2 . However, 
the Figures seem to indicate that any scaling of l± between these two extremes leads to a 
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leading order contribution and has a pinched pole. The question then arises why we do not 
have to consider more modes, such as one with momentum scaling as ~ (A 2 , A 2 , A 3 / 2 ). 
Expanding around this point, we would obtain a mode in the effective theory similar to the 
construction discussed above. The important point, however, is that the overlap contributions 
between this mode and the Glauber and soft modes are identical to the new mode itself. Thus, 
after taking the overlap contributions into account, the contributions from any additional 
modes vanish. A similar argument can be made that no additional modes are required even 
if power suppressed terms are considered. 

3.2 Relevance of Glauber gluons to other processes 

In this Section we will assume that a generic one-loop diagram has a pinch singularity in 
the Glauber region and derive some necessary characteristics of the diagram from this and 
a few other assumptions. The following analysis is closely related to the preceding Section. 
Previously we started by studying a loop diagram with given external leg momenta and found 
that the Glauber region is pinched. Here we are going the opposite way: Let us start with 
the assumptions that 

1. the loop momentum I ~ (A 2 , A 2 , A) leads to a pinch singularity, 

2. external momenta hi are on the mass-shell, 

3. a momentum component can be at most 0(1), no inverse powers of A. 

Our choice will be to perform the l + integration by contours. Consider the propagator 
(I + k)~ 2 , which leads to a simple pole of the form (l + — z)~ l with 

,-<k+fa£z£-»+. (3.5) 

We will use the following ansatz for the scaling of k, namely /c ± ~ 0(A ra± ), k± ~ 0(A ra± ). 
The on-shell condition means that n+ + n_ = 2n±. Since / + ~ 0(A 2 ), it follows that z must 
also scale like 0(A 2 ), and one of the following conditions (a) or (6) must be satisfied: 

(a) 2min(l,n_i_) — min(2, n_) = 2 and n + > 2 , 
(6) 2min(l,nj_) — min(2, n_) > 2 and n + = 2 . 

It is easy to show that condition (6) can never be satisfied with the assumption that n± , n± > 
0, and that the only scaling of k consistent with (a) is k ~ (A 2 , 1, A), i.e. collinear. Strictly 
speaking we get n_ = and n± > 1, e.g. (A 4 , 1, A 2 ) would also be allowed. The momentum 
with largest invariant mass, however, (the least restrictive condition) is the standard collinear 
one. External legs with even smaller invariant masses will pinch the loop integral also in a 
Glauber- type configuration, for example / ~ (A 4 , A 4 , A 2 ), which is readily rescaled by A' = A 2 
to yield the same result. 
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Figure 4: A one-loop di- 
agram with loop momen- 
tum I. 



We have therefore found that attaching an incoming collinear 
leg to the loop is consistent with a Glauber loop momentum, al- 
though certainly not sufficient. In particular, note that k~ domi- 
nates over l~ . In order to obtain a pinch singularity, we will hence 
need another leg with collinear scaling, but with outgoing momen- 
tum 5 . Next, we note that it was a choice to perform the l + integra- 
tion by contour, which pins down the / + scaling but leaves /~ open. 
We therefore repeat the above argument starting with the l~ com- 
ponent by contours and find that we need to attach two more legs 
to the loop diagram, one in- and one outgoing with anti-collinear 
momentum scaling. 

It follows that any process involving diagrams of the form depicted in Figure || will 
be sensitive to Glauber interactions. The characteristics are that there are at least four the 
external legs, with two of them collinear and two of them anticollinear. For a pinch singularity, 
each pair must consist of an incoming and an outgoing momentum. The simplest such case is 
forward scattering of a collinear with an anti-collinear field, where no other external legs are 
attached to the diagram. Multiple exchanges of Glauber gluons will give rise to an effective 
potential between the collinear and anti-collinear fields when integrating Glauber gluons out 
of the theory. This is similar to the exchange of potential gluons in non-relativistic QCD, 
and can be found elsewhere in the literature [^9|, |50| . In our case of Drell-Yan far from 
machine threshold there is one extra leg in the diagram from which the lepton pair will 
emerge. From the point of view of Glauber gluon relevance this case is identical to Higgs 
production, or indeed any New Physics particle production in hadronic collisions away from 
machine threshold. 



4. Conclusions and Outlook 

Using an explicit matching calculation we have investigated whether Glauber gluons are re- 
quired as a degree of freedom in SCET. Focusing on the operator O2, whose Wilson coefficient 
is well known, we showed that Glauber gluons are required if the matching is performed with 
final states that contain both initial and final collinear particles in the same collinear direction. 

Even though we did not directly consider the Drell-Yan process, our choice of external 
states was such that all contributions of the Drell-Yan process, namely spectator-active and 
spectator-spectator interactions in addition to active-active were present. This allows us to 
make parallels between our conclusion about the consistency of matching and the correct 
modes for the Drell-Yan amplitude. Our conclusion is that for the exclusive Drell-Yan am- 
plitude the correct effective theory would require Glauber modes. Note that we have not 
discussed under what circumstances the contribution of Glauber gluons cancel when squaring 
the amplitude. 

For our analysis it was important to avoid double counting between the modes by per- 

forming zero-bin subtractions f48|j fro m the collinear and Glauber modes. It is worthwhile to 

5 The next propagator will then be (I + K)~ 2 with K = k — k! . Since both k and k' have collinear scaling, 
so does K. 
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emphasize the spectator-active case for SCETq, since in this topology there is an interesting 
example of non-vanishing zero-bin subtraction from the naive collinear mode (-Z4 ° )og 7^ 0. 
It shows that taking these subtraction serious renders the effective theory robust toward the 
introduction of other modes and does not invalidate previous calculations. 

Prom the scaling of the Glauber mode it is apparent that it cannot be made exactly on- 
shell. However in the case of the spectator-spectator topology there is a pinch singularity in 
the Glauber region. Thus we find an apparent contradiction to the Coleman-Norton theorem 
[O], which states that all pinched surfaces arise from on-shell degrees of freedom. This 
contradiction is resolved if one notices that Landau Equations used in the Coleman-Norton 
theorem are conditions for appearance of a true singularity, which correspond to the case 
when the loop integral is infinite, i.e when power counting parameter of the effective theory 
vanishes: A = 0. In this limit both the Glauber and soft momenta are identical to each other 
Zg = Is = 0. Of course in this limit both modes become on-shell, so the contradiction with 
Coleman-Norton theorem goes away. 

It would be interesting to reformulate the Landau equations taking into account the 
proper power counting, i.e. instead of writing down a condition of having a true pinch 
singularity, find a condition for pinched poles to occur at distance of say order A 2 from 
each other. It should be possible to use such an analysis to discover both Glauber and 
soft pinches in the spectator-spectator diagram directly from such relaxed Landau Equations. 
Understanding this question in details might potentially be of practical importance since it can 
lead to a nice recipe for arbitrary process on how to read off the correct long distance modes 
(valid to all orders in perturbation theory with a possibility to include power corrections). 
This, however, is beyond the scope of the present paper and will be studied elsewhere. 

Another important step is to study the effect of Glauber gluons in SCET on physical 
observables, and not just amplitudes as done in this paper. The expectation here would 
be to understand the cancellation of Glauber gluons in the inclusive in the final hadronic 



state Drell-Yan cross-section, which was proved in full QCD in Refs. |16|, 17, 18|. The main 
challenge in doing so is that the Glauber mode scaling is such that the corresponding particle 
is always off-shell, such that it can not be naively included in the SCET Lagrangian. A better 
way to proceed might be to interpret this mode as an effective potential, giving rise to forward 
scattering of two collinear particles in different directions [46|. 
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